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ABSTRACT 

We have previously analyzed sensitive mid-infrared observations to establish that the Pipe 
Nebula has a very low star-formation efhciency. That study focused on YSOs with excess infrared 
emission (i.e, protostars and pre- main sequence stars with disks), however, and could have missed 
a population of more evolved pre-main sequence stars or Class III objects (i.e., young stars 
with dissipated disks that no longer show excess infrared emission). Evolved pre-main sequence 
stars are X-ray bright, so we have used ROSAT All-Sky Survey data to search for diskless pre- 
main sequence stars throughout the Pipe Nebula. We have also analyzed archival XMM- iVewton 
observations of three prominent areas within the Pipe; Barnard 59, containing a known cluster of 
young stellar objects; Barnard 68, a dense core that has yet to form stars; and the Pipe molecular 
ring, a high-extinction region in the bowl of the Pipe. We additionally characterize the X-ray 
properties of YSOs in Barnard 59. The ROSAT and XMM- iVewton data provide no indication 
of a significant population of more evolved pre-main sequence stars within the Pipe, reinforcing 
our previous measurement of the Pipe's very low star formation efficiency. 

Subject headings: stars: pre-main sequence — X-rays: stars 



1. Introduction 

In a study using deep mid-infrared imaging, we 
found a very low level of star form ation activity in 
the P ipe Nebula (10'* Mq, 130 pc. lLombardi et al 



20061 ) ■ with an upper limit of only 21 candidat e 



YSOs in the entire cloud ()Forbrich et al.l |2Q09|) 



This survey was most sensitive to young stellar 
objects (YSOs) with infrared excesses and presum- 
ably in the earliest stages of evolution. A popu- 
lation of more evolved class III YSOs, i.e., disk- 
less young stellar objects without infrared excess 
(e.g., iLada 1987 ). if present, would have impor- 
tant implications for the cloud. Most importantly, 
such a population would lead to an increased esti- 
mate of the region's star-formation efficiency and 
would also indicate that the molecular cloud is 
more long-lived than otherwise thought. 

Class III objects can be identified by their X-ray 
emission. While this technique has been known for 



some time (e.g.. IKu fc Chaiian Il979l) . it has been 
strengthened by the analysis of the Chandra Orion 
Ultra-deep Project which produced a statistically 
well de fined description of the X-ray properties of 
YSOs (jPreibisch fc Feigelso n 2005). Furthermore, 
class III objects have been studied in Spitzer in- 
frared surveys of othe r nearby (d ^ 1 20-16 pc) 



star-forming regions. iPadgett et al. (|2006l) tar- 



geted 83 class III sources in the Taurus, Lupus, 
Ophiuchus, and Chamaeleon star-forming regions 
at 24 /xm, detecting all 83 but identifying excesses 
in only 5. ICieza et"aL ( 2007 ) also used c2d data 
to study 102 spectroscopically identified class III 
sources in the Lupus and Ophiuchus regiontQ, de- 
tecting 80% of them at 24 (82/102) at a sensi- 
tivity comparable to our data and with distances 
of about that of the Pipe. Requiring mid-infrared 
(MIPS) as well as X-ray detections thus provides 
better constraints on class III source populations 
than any of these two datasets alone could pro- 



^now at: Department of Astronomy, Cornell University, 
Ithaca, NY 14853, USA 



^leaving out the more distant IC 348 in Perseus for the sake 
of this comparison 
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vide, and we use this additional constraint in our 
analysis of the Pipe Nebula. 

In this study, we first use the ROSAT All-Sky 
Survey (RASS) as the only available means to 
search the entire Pipe Nebula complex for X-ray- 
active stars. We then study three regions of partic- 
ular interest using higher-sensitivity data obtained 
by the X-ray satellite XMM-Newton. The latter 
satellite has considerably more collecting area and 
is sensitive to harder X-rays than ROSAT, allow- 
ing observations of more extincted sources. We use 
XMM-JVewton to probe three prominent regions in 
the Pipe Nebula complex in more detail; Barnard 
68 (B 68), Barnard 59 (B 59), and the Pipe molec- 



Spitzer-MlPS study (jPorbrich et al.ll2009l) 



ular ring in the bowl of the Pipe (Muenc h et al 



20071) . in the following referred to as the PMR. 
While B 59 is known t o contain a clust er of YSOs 
(|Brooke et al., 2007: Covev et al.ll201o( ). the entire 
bowl area has been found to be devoid of proto- 
stars ( Forbrich et al. l2009l) . even though the av- 
erage extinction in that area is at levels of visual 
extinction higher than those of B 59. 

The Pipe region is seen in projection against 
the Galactic bulge, providing an enormous back- 
ground source density, particularly at near-infrared 
wavelengths. Therefore, we have to compare 
the density of X-ray sources inside and outside 
of the identified den se core population (from 



Rathborne et al. 20091 ) to estimate the local den- 
sity of background X-ray sources. In the follow- 
ing, we use the same extinction criterion in defin- 
ing the 'on-core ' region s that was also used by 



Rathborne et al. (I2n09h : These regions have a 
minimum optical extinction, Av{snh)> 1.2 mag, 
in an extinction map in whic h a slowly varyiii g 
background has been removed (|Alves et al.ll2007t) . 
We ha ve used the same cr iterion in our Spitzer 
study (jPorbrich et al.ll2009l ). For purposes other 



than defining the on-core regions , we use the to- 



tal ex tinction map computed by iLombardi et al 
(|2006l ). The pixel size of the extinction maps is 
30". 

Looking for possible class III objects, we se- 
lect X-ray sources with 2MASS counterparts in 
these regions, requiring Ks band detections. We 
then use their X-ray properties as well as Spitzer 
IRAC and MIPS data to better characterize them. 
Finally, to further constrain any population of 
class III sources, we combine the X-ray source se- 
lection with detections at 24 /im from our sensitive 



Observations 



To cover the entire Pipe Nebula complex, 
we make use of the ROSAT All-Sky Survey 
in the form of the Bright Source Catalogue 



(IVoees et al 



of t he 
I999I1 



, and the Faint Source Cata- 
logue (jVog es et alllioooh . The observations for 



this survey were carried out in 1990/91 with the 
Position- Sensitive Proportional Counter (PSPC), 
in the energy range of 0.1-2.4 keV. 

The X-ray satellite XMM-JVewton has observed 
the Pipe Nebula region several times. Here, we an- 
alyze three particular observations, targeting the 
prominent extinction cores B59, FeSt 1-457g, and 
B68. The latter region was observed for 51 ksec 
on September 27, 2002 (07:19-21:00 UT, median 
mter, observation ID 0152750101), followed by 
observations of the two other regions on August 
29, 2004: Fest 1-457 (06:19-19:53 UT, 53 ksec, 
thin fiher, observation ID 0206610201) and B59 
(21:53-10:14(-H1) UT, 49 ksec, thin filter, obser- 
vation ID 0206610101). Ah of these observations 
hav e been processed as part of the 2XMMi cata- 
log ( Watson et al.ll2009l) . Note that pipeline prod- 
ucts like spectra and light curves for all catalog 
sources can be downloaded from a dedicated web- 
sitdl. The good time intervals for these obser- 
vations amount to observing times of 48.2 ksec 
for B59, 43.8 ksec for the PMR, and 44.8 ksec 
for B68. The energy range accessible to XMM- 
Newton and used for the 2XMMi catalog is wider 
than the one of ROSAT; it spans 0.2-12 keV. In 
the following, we mainly work with this catalog, 
but we have also reprocessed the data using SAS 
8.0.0 to subsequently run our own source detec- 
tion using the SAS task ewavelet. We also employ 
the SAS task especget to produce X-ray spectra 
of selected sources. For this purpose, we mostly 
use the EPIC-pn data only, rather than EPIC- 
MOS data, apa rt from the unres olved sources 
[BHB2007] 1&2 (jBrooke et al.ll2007t) which are so 
far off-axis {r = 11. '9) that they are only detected 
with the EPIC-MOS detectors. Typically, spectra 
were extracted from areas with a radius of 25" 
and similarly sized, nearby background regions. 



^Note that although this observation is hsted as targeting 
FeSt 1—457, it is pointed at the neighboring PMR. 



^e.g. , ilhttp : //xcatdb . u- strasbg ■ f r/2xnmi | 
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The only exception is the close pair [BHB 2007] 
6 & 7, where we chose a radius of 10" in an at- 
tempt to separate the two components. The task 
especget produces the corresponding ancillary re- 
sponse and response matrix files. Spectra were fi- 
nally binned with minimum bin sizes of 20 counts. 
Fits were carried out using the CIAO tool Sherpa, 



V3.4 ( Freeman et al.ll2001^ 



In order to assess the infrared spectral en- 
ergy distributions of selected sources, we have 
performed aperture photometry on Spitzer-IRAC 
data covering the region, most notably observation 
IDs 5132288, 14089472, and 14090496. For this 
purpose, we have used pipeline PBCD data from 
the Spitzer archive. The typical aperture radius 
was 3 image pixels with an annulus for background 
subtraction with inner and outer radii of 10 and 20 
pixels. We have ad ditionally used Spitee r-MIPS 
data as discussed in Forbrich et al. I (l2009l) . 



3. Analysis 

3.1. The ROSAT All-Sky Survey 

The area covered in our mid-infrared survey 
(|Forbrich et al.ll2009l Fig. 8) contains 44 ROSAT 
All-Sky Survey sources, i.e., 12 and 32 sources 
from the Bright and Faint Source Catalogues, re- 
spectively. The locations of these sources are 
shown in Figure [TJ overlaid on an extinction map 
of the Pipe Nebula. Given typical la positional 
uncertainties of ~ 20", it is nearly impossible to 
reliably identify near-infrared or optical counter- 
parts to ROSAT X-ray sources in this region with 
its high density of background sources, simply due 
to confusion. While this is somewhat easier at 
mid-infrared wavelengths due to the lower source 
density in this wavelength range, all 36 sources 
whose 3cr error circles are located entirely within 
the MIPS source detection coverage have one or 
more mid-infrared counterparts at 24 /im. Eight 
additional sources have error circles that only 
partly overlap with the area where MIPS source 
detection was performed. Note that the error cir- 
cle of only one ROSAT source (IRXS J173319.3- 
255416) overlaps with the on-core regions and all 
but two of the sources have zero extinction in the 
extinction map where the extended structure has 
been removed (see above) . Their median total ex- 
tinction is Ay = 2.4 mag. Correlating the source 
list with catalogs using SIMBAD shows that two 



of these ROSAT sources have been identified as 
X-ray binaries, two as Hipparcos stars (at 54 and 
149 pc, respectively) and one as an RS CVn bi- 
nary. 

To assess the completeness of the ROSAT All- 
Sky Survey for our purposes, we determine the 
limiting unabsorbed X-ray luminosities for differ- 
ent levels of foreground extinctions. The limit of 
six source counts in the Faint Source Catalog cor- 
responds to a limiting count rate of 0.02 s^^ with 
the typical exposure time of 300 s in this area. 
We can now simulate X-ray sources with a given 
spectrum to find the limiting unabsorbed X-ray 
fluxes and in turn the luminosities corresponding 
to this count rate limit. Using the Portable . Inter- 
active Multi-Mission Simulator (PIMMfl iMukail 



1993|), we determine such unabsorbed X-ray lu- 



minosities for sources using simulated Raymond- 
Smith thermal plasma spectra with temperatures 
of ~ 10 MK. While the median total visual extinc- 
tion from our extinction map for the 44 ROSAT 
sources is Ay =2.4 mag, we calculate the limit- 
ing luminosities also for Ay =0,5, and 25 mag, 
the maximum of our total extinction map (which 
constitutes only a lower limit for the actual max- 
imum extinction in the cluster). The visual ex- 
tinctions can be converted into absorbing hydro- 
gen column densities using the empir ical relation 



iVff [cm-^jRj 2 X 10 ^^ X Ay [mag] (e.g. 



Rvteilll996l: 



Vuong et al. 2003[ ). Conservativ ely assuming a 



plasm a temperature of 1 keV fe.g.. lFeigelson et al 



20051 ). the limiting luminosities are \og{Lx)=29.5, 
30.2, 30.7, and 32.3 erg s^^ for the given visual ex- 
tinctions (in rising order), assuming a distance of 
d = 130 pc. Note that the high-extinction regions 
cover only a minute fraction of the complex. 

We can now make use of the cumulative X- 
ray luminosity distribution functions derived for 
the PMS population of the O rion Nebula Cluster 
( Preibisch fc Feigelson 20051 ). to assess the com- 
pleteness of the ROSAT data for G, K, and M pre- 
main-sequence stars. The limiting unabsorbed X- 
ray luminosity that we have determined for the 
median foreground extinction of Ay = 2.4 mag 
corresponds to a completeness of ^ 55% for G and 
K stars and ^ 15% for M stars. For Ay — 5 mag, 
the completeness diminishes to ^ 30% for G and 
K stars and close to zero for M stars while with- 



• ^http: //heasarc .nasa. gov/docs/sof tware/tools/pimms .htmll 
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out foreground extinction, the numbers are ~ 90% 
completeness for G and K stars and ~ 50% for M 
stars. 

A significant subset of tire selected ROSAT 
sources can be expected to be located in the Galac- 
tic bulge against which the Pipe Nebula is pro- 
jected, to be of extragalactic nature, or to be in the 
foreground. To assess this aspect, we determine 
the source density in the Pipe and in neighbor- 
ing regions (shown in Fig. [T|). The overall RASS 
source density in the region covered in our mid- 
infrared survey, where we find 44 RASS sources, 
is 3.29 sources per square degree. However, it 
turns out that the RASS source density in the 
Pipe Nebula is lower than in the immediately sur- 
rounding regions. Our Pipe Nebula total extinc- 
tion map covers nearly 49 square degrees which is 
about 3.7 times more than the area covered in the 



mid-i nfrared survey (see Figure 8 in lForbrich et al 
2009t ). This area contains a total of 199 RASS 
sources. In the immediate surroundings of the 
high-extinction regions inside the infrared survey 
area, we thus find a RASS source density of 4.37 
sources per square degree. This source density is 
a good measure of the contaminant source density 
in this region since the median total extinction of 
the RASS sources inside the infrared survey area 
[Ay = 2.4 mag) and the median extinction for 
RASS sources in the remainder of the area covered 
by the extinction map {Ay = 2.0 mag) are very 
similar. Extrapolating from the source density in 
the surroundings, we would expect to find 58.4 
sources in the infrared survey area, but we only 
find 44 of which six are known to be unrelated to 
the Pipe Nebula. The Poisson error ranges of these 
two numbers do not overlap. It appears likely that 
the source density in the Pipe Nebula is lowered 
by the partial extinction of a distribution of back- 
ground sources, for example in the Galactic bulge 
or extragalactic sources, by the molecular cloud. 
This effect is difficult to quantify since it depends 
on the intrinsic source spectra, e.g., on whether a 
given class of sources is rather soft and therefore 
more susceptible to extinction. However, while in 
our total-extinction map of the Pipe Nebula and 
its surroundings, 11.5% of the pixels are of extinc- 
tions Ay > 5 mag, only 6.8% of the RASS sources 
in the same area are affected by the same amount 
of extinction, as derived from the same extinction 
map, consistent with a picture in which we do not 



detect the most heavily extincted sources. 

Without sufficient information about optical or 
infrared counterparts, we can only make use of 
the X-ray source properties for further characteri- 
zation. For the ROSAT sources, leaving aside the 
potentially variable count rates, these are primar- 
ily th e two standardized hardness ratios HRl and 
HR2 ( Voges et al. 19991 ). ranging from -1 to -|-1. 
Similar to, e.g., infrared colors, the hardness ra- 
tios provide a crude characterization of spectral 
properties in the absence of actual spectra. The 
second ratio probes harder X-ray emission within 
the ROSAT band than the first hardness ratio. 
However, among the 44 ROSAT sources selected 
here only six have hardness ratios HRl and HR2 
with errors of < 0.3. 

In order to judge even these few hardness ra- 
tios, we need a reference region. Since the Taurus 
star-forming region is of about the same size, mass 
and distance as the Pipe Nebula, we compare our 
sample with a correlation of Taurus PMS stars and 
the ROSAT All-Sky Survey, using the listing of a ll 
known Taurus PMS stars bv lKenvon et al.l (|2008l ). 

While this list is not complete, it provides us 
with a representative comparison sample. Con- 
servatively using 2(7 error circle^, there are 46 
RASS sources encompassing at least one PMS ob- 
ject from this Taurus list. This overall number 
of sources is simi lar to the results of the ROSAT 
Taurus survey bv lNeuhauser et aL ( 19951 ). Out of 
these 46 sources, 29 have hardness ratios with rea- 
sonable errors < 0.3. In Fig. [2l we show a plot of 
the two hardness ratios for the Pipe and the Tau- 
rus samples. While the latter sou rces delineate the 



same region of the diagram that iNeuhauser et al 



()l995l) bund to be occupied by class II and HI 
objects, only about half of the Pipe sources are in 
that region, although the overall numbers are very 
low. Unfortunately, PMS objects are not the only 
source class with predominantly positive hardness 
ratios: X-ray binaries fall onto the same area of 
the hardness ratio diagram (e.g.. lHaberl fc Pietsch 
1999t iHaberl et ailbOOOl) . To show only one ex- 



ample: While both the only ROSAT source in the 
on-core region and the brightest ROSAT source in 
our sample have HR1=HR2=1, the latter is a pre- 



^ Using Icr errors would yield 24 sources, while a 3cr search 
radius yields 56 sources. The median 2ct search radius is 
28". 
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viously known X-ray binary that is unrelated to 
the Pipe Nebula. As a further complication, the 
intrinsic hardness ratios of background sources are 
affected by foreground extinction. The hardness 
ratios therefore do not provide significant further 
constraints for our purposes. 

In summary, the ROSAT source density in our 
Pipe Nebula mid-infrared survey area is lower than 
in the immediate surroundings; 39 of 44 sources 
remain unidentified. The on-core regions con- 
tain at most a single RASS X-ray source. Of 
all 44 sources, their median total extinction is 
Av =2.4 mag, only six have hardness ratios of suf- 
ficient quality to compare them to Taurus PMS 
objects, and only about half of them fall onto the 
same region of the hardness diagram as the Tau- 
rus sources, even though that parameter space of 
hardness ratios is not unique to PMS sources. We 
conclude that likely none of the RASS sources lo- 
cated toward the Pipe Nebula are class III YSOs. 
To substantiate this result, we take a deeper X-ray 
look at three high-extinction regions in the Pipe 
Nebula using XMM-JVewton. 

3.2. XMM-Newton data 

As noted above, all of the XMM-iVewton ob- 
servations considered here have been_processcd as 
part of the 2XMMi catalog fjWatson et al. 2009). 
We use this catalog to search for X-ray sources 
with 2MASS- JC5 band counterparts, particularly 
in the on-core regions. Ancillary Spitzer MIPS 
and IRAC data allow us to subsequently study the 
SEDs of identified sources in more detail. Only 
one ROSAT source is located within the fields 
of view of the three XMM-JVewton observations 
(IRXS J172236.0-235938 near the southern edge 
of the B 68 field; it is compatible with the position 
of a 2XMMi source). 

3.2.1. Initial X-ray source selection 

In a first step, we exclude low-quality detections 
from the X-ray sources listed in the 2XMMi cat- 
alog for the three different observations. After a 
visual check, we exclude five sources listed as likely 
spurioufH. Out of 297 initially listed sources, this 
leaves 292. Only a minor filtering of detection like- 
lihoods has been applied; we conservatively use a 



^We exclude sources with a 2XMMi summary flag of 4. 



minimum detection likelihood of srcML=8, as pro- 
vided by the 2XMMi catalog (where the minimum 
detection likelihood reported is srcML=6). We ex- 
clude 17 sources that were detected as extended 
objects, identified by their catalogued extension 
likelihood. These two steps of filtering leave 251 
X-ray sources. The distribution among the three 
observations is as follows: 112 of these sources are 
located in the field of view of the B 68 observation, 
101 in the field of B 59, and 38 in the PMR. 

Any list of X-ray detections will likely include 
background sources in the Galactic bulge, and 
we note that the three fields covered by XMM- 
Newton cover a total of 43704 2MASS sources. We 
therefore again separate parts of the analysis into 
on-core and off-core regions. It is important to 
keep in mind that the on-core regions cover dif- 
ferent fractions of the three fields of view. While 
these fractions are 49% and 39% for B 59 and the 
PMR, respectively, only 14% of the B 68 field of 
view consists of on-core regions. The average total 
extinction in the three fields of view, as measured 
in our total extinction map, is highest in the PMR 
{Av = 8.7 ± 2.6 mag), followed by a similarly ex- 
tincted B 59 [Ay — 5.4 ± 3.5 mag) and a quite 
different B 68 area at only Ay = 1.2 ± 1.2 mag. In 
contrast to the other two regions, there is virtually 
no extinction surrounding the B 68 core. These 
basic properties of the three fields as well as re- 
sults of the subsequent analysis are summarized 
in Table m 

Among the X-ray properties, we only make use 
of the X-ray count rates in the following. The 
significant extinction influences the intrinsic hard- 
ness ratios of background sources so that hardness 
ratios from reference regions cannot be directly ap- 
plied. As we have seen in the discussion of ROSAT 
data, it is nearly impossible to identify PMS ob- 
jects based on their X-ray hardness ratios alone. 
Therefore, the X-ray hardness ratios provided by 
the 2XMMi catalog do not allow us to constrain 
source classes. 

3.2.2. Infrared counterparts and SED analysis 

In a second step in our search for candidate 
class III sources, we only select X-ray sources 
with 2MASS counterparts and, more specifically, 
i^s-band detections. To define a cutoff Ks 
band magnitude, we use the fact that 90% of 
the unique 2MASS counterparts to the Taurus 
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PMS sources listed by iKenvon et al. ( 2008 ) have 
Ks < 12.4 mag, using all cases where the source 
is within the 3a 2MASS error circle. While the 
source density in the Pipe Nebula region is very 
high due to its projection against the Galactic 
bulge, we can thereby at least eliminate some X- 
ray sources without sufficiently bright NIR coun- 
terparts. Note that the total number of 2MASS 
sources fulfilling this Ks band cutoff criterion is 
still 7630 for the three pointings. 

For every 2XMMi source selected thus far, we 
have extracted the X-ray positions corresponding 
to the three observations that we consider and we 
look for 2MASS counterparts in the Scr error cir- 
cles as quoted by the 2XMMi catalog. The ra- 
dius of the median 3a error circle is 3'.'9 where 
the values range between 1'.'6 and 8'.'4. Start- 
ing with the list of 251 X-ray sources, we thus 
find 46 X-ray sources with 2MASS counterparts 
of Ks < 12.4 mag. One of the sources has two 
suitable 2MASS counterparts in the search radius 
in which case we choose the one closer to its X-ray 
position. Note that the positions of 10 out of these 
46 X-ray sources are also compatible with weaker 
2MASS sources. 

Given the high 2MASS source density in the 
target regions, we can now assess the possibility 
of chance alignments, i.e., the number of 2MASS 
sources above our Ks band cutoff that we would 
find if both the X-ray search areas and the 2MASS 
sources were randomly distributed across the fields 
of view. A summarized in Table [1] we would glob- 
ally expect about 17 of the 46 selected sources to 
be chance alignments. Not differentiating between 
on- and off-core regions, only the number of se- 
lected sources in the B 59 and B 68 fields is signif- 
icantly above the expected number of such chance 
alignments. 

To find out whether any of the 46 selected X-ray 
sources with near-infrared counterparts, includ- 
ing the subgroup of 13 sources in the on-core re- 



^The 2XMMi catalog uses a systematic error of either 0'.'35 
or I'.'OO, depending on whether an acceptable astromet- 
ric correction was obtained. Among the datasets consid- 
ered here, the smaller error applies only to the B 59 ob- 
servation. The infrared counterparts of two known YSOs 
in B 59, imaged as bright sources at large off-axis angles 
(r = 11. '9) and including one potential unresolved double 
X-ray source, lie just outside the total 3cr error circles as 
listed in the 2XMMi catalog. We therefore assume a sys- 
tematic error of 0'.'50 for B 59 instead of 0'.'35. 



gions, have spectral energy distributions indicative 
of their YSO status, we have obtained Spitzer mid- 
infrared aperture photometry for these sources, us- 
ing the four IRAC bands in addition to our pre- 
viously discussed MIPS 24 /im data. While any 
candidate class III sources might simply show pho- 
tospheric SEDs, we here primarily look for previ- 
ously missed excess sources. All selected X-ray 
sources in on-core regions and most in the off-core 
regions were covered in all IRAC bands and at 
24 /xm and 70 /xm (MIPS). Apart from the known 
YSOs, all X-ray sources with NIR counterparts 
and Spitzer photometry were su bsequently fitted 
using the online SED fitter by iRobitaille et al 



(|2007l ) to check both YSO models and extincted 
stellar photosphere fits. The model fits reveal that 
all but one of these sources have SEDs that can 
be reasonably explained by stellar photosphere^. 
The extinction derived from these SED fits is, in 
most cases, lower than the estimate based on the 
extinction map by 1-2 mag, indicating that the 
objects are located inside the cloud. In the case of 
source 5 (2XMMi J171117.5-272535), the SED fit 
requires an extinction of only Ay = 0.3 mag while 
the corresponding pixel of the extinction map in- 
dicates Av = 13.0 mag. A fit to the X-ray spec- 
trum supports a low absorbing column densitjl^, 
as do the soft hardness ratios of this source, as 
catalogued. This source, therefore, probably is a 
foreground star. 

3.2.3. Final source selection and assessment 

To better constrain potential candidate YSOs, 
we can apply a lower limit to the X-ray fluxes, 
knowing that X-ray counterparts to any YSOs 
would be relatively bright X-ray sources, at 
least in extinction-free regions. We thus ex- 
clude all sources with unabsorbed X-ray fluxes 
that are lower than those derived from the X- 
ray lu minosities of 90% of the O rion PMS M- 
stars (jPreibisch fc FeigelsonI 120051) . i.e.. Fx < 
2.5 X 10"^'' ergs~^cm~^. In order to determine 



*The exception is 2XMMi J173400.2-253746 where only the 
24 ^m flux is significantly brighter than photospheric (the 
source is undetected at 70 fim). This source is outside of, 
but close to the PMR on-core region. 

9 While an absorbed APEC fit to 101 total EPIC-pn counts 
does not produce statistically meaningful parameter es- 
timates, it yields a 3a upper limit for the extinction of 
Av < 0.7 mag, using the above-mentioned conversion. 
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the unabsorbed X-ray flux for a given source, we 
use PIMMS to convert the total count rate (0.2- 
12keV), conservatively assuming that the spec- 
trum is thermal plasma emission that can be ex- 
plained by a 1 keV APEC model and the extinc- 
tion from the total extinction map. Using this 
criterion, we can exclude nine low-extinction (and 
off-core) objects of the 14 sources in the B 68 field, 
reducing the sample to a total of 37 X-ray sources 
with 2MASS counterparts. The X-ray sources in 
the surroundings of the B68 core appear to be 
unrelated to the cloud. 

Of these 37 sources, 23 sources are located in 
the B 59 field, 9 in the PMR, and 5 in the B 68 
field. A view of the three X-ray datasets and the 
selected sources is shown in Figures [Sj lU and [5] 
Among these are 13 sources in the on-core regions 
of B59 (12) and the PMR (1), with none in the 
on-core regions of the B 68 field. These sources 
are listed in Table [2J The single source in the 
PMR and one source in B 59 are close to the Ks- 
band cutoff. Only the number of on-core sources 
in the B 59 fields surpasses the number of expected 
chance alignments, as summarized in Table[T] The 
difference between B 59 and the PMR is striking, 
given that both are comparable in overall extinc- 
tion and amount of on-core regions covered. 

Among the 13 X-ray sources with 2MASS- 
Ks counterparts in the on-core regions of B 59, 
there are 7 prev iously identified YSOs from 
Brooke et all (l2007t ). not differentiating [BHB2007] 



1 and 2 which are unresolved in the XMM- JVewton 
data due to the large off-axis angle. Additionally, 
[BHB2007] 7 appears to have a weak X-ray coun- 
terpart that is not listed in the 2XMMi catalog, 
but was identified in our wavelet source detection. 
One X-ray detected YSO - [BHB2007] 20 - is lo- 
cated outside of the on-core region, and the deeply 
embedded protostar [BHB2007] 11, while X-ray- 
detected, is fihered ou10 due to its weak 2MASS 
counterpart. In the innermost, most deeply em- 
bedded part of the S pitzer-identified young cluster 



( Brooke et al.l 120071 ) . no new candidate members 



were detected. However, a 2MASS source south 
of [BHB2007] 6 and 7 (not among the identified 
YSOs) was found to have a weak, soft X-ray coun- 
terpart. The X-ray properties of the known YSOs 



in B 59 are discussed in more detail in Section [3.3l 
Since we undoubtedly again detect background 
sources in the Galactic bulge next to putative 
sources in the cloud as well as foreground sources, 
we again analyze the density of X-ray sources in 
the on-core and the off-core regions. In B 68, all 
X-ray sources with NIR counterparts lie in off- 
core regions; this is not the case in B 59 and the 
PMR. In the absence of a class III source popula- 
tion in the Pipe region, we would expect to find a 
source density that is lower than or equal to the 
source density in the off-core region due to the ad- 
ditional extinction. In B59, with its known clus- 
ter of YSOs, we find that the source density in the 
on-core regions is indeed slightly higher than the 
source density in the surrounding off-core regions 
(52% of the X-ray sources with NIR counterparts 
lie in the on-core regions which cover 49% of the 
field of view). In the PMR, in striking contrast, 
we find that only 11% of these sources are in the 
on-core regions even though these cover 39% of the 
field of view. Note that the source density in B 59 
turns into an underdensity without the counter- 
parts of known YSOs. This comparison of X-ray 
source densities in the on-core and off-core regions 
shows that there is no indication for a significant, 
previously unknown population of PMS objects. 

In the on-core regions of the XMM- JVewtoii ob- 
servations, we therefore find an upper limit of 6 
new evolved YSOs since this is the number of X- 
ray sources with appropriate -ftTs-band counter- 
parts, not counting the known YSOs. However, 
the fact that the two regions with any detections 
in the on-core regions, B 59 and the PMR, both 
show underdensities with respect to the surround- 
ing background source densities (not taking into 
account the known YSOs in B 59) suggests that a 
more realistic upper limit is closer to zero. A more 
stringent upper limit results when additionally re- 
quiring detections at 24 /im, assumi ng that this 
selects 80 % of the class HI sources (Cieza et al 



^"Note that in our search for class III sources, the criteria 
were not designed to include protostars. 



i2007L see above). Apart from the known YSOs, 
only two X-ray sources with Ks-hai\A detections 
in the on-core regions have 24 /im counterparts, 
both of them in B59, suggesting an upper limit 
of only a handful of class HI sources. The com- 
bination of the three selection criteria (2MASS, 
X-ray fiux, and MIPS detection) selects at least 
two thirds of any class HI population. 

Finally, we can additionally give strict upper 
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limits for any evolved YSO population in the 
three entire fields to account for any more dis- 
tributed populations. For this purpose, we cannot 
treat the off-core source density as entirely un- 
related to the cloud. While we have found that 
all sources in the B 68 and ring fields can be ex- 
plained by a combination of chance alignments and 
weak background sources, there are 23 selected X- 
ray sources in the B 59 area when only 3.5 would 
be expected from chance alignments. Subtract- 
ing the expected number of chance alignments and 
the number of known YSOs among the selected 
sources (7), we derive an upper limit of ^ 13 new 
candidate YSOs in the entire B 59 XMM-JVewton 
field. Note that among the selected sources, the 
field contains six detections at 24 /im apart from 
the known YSOs, again suggesting a more modest 
upper limit for any additional more evolved YSOs. 
These sources would be in addition to the upper 
limit of 15 previously known candidate YSO s in 
B59 (jBrooke et al.ll2007t |pb7brich et al.ll2009[) . 



3.3. Barnard 59 

In addition to the analysis in the previous sec- 
tion, which focused on finding candidate YSOs 
other than the ones already known from mid- 
infrared imaging, we now take a more detailed 
X-ray look a t the s ample of YSOs identified by 
Brooke et al.l (|2007l) . using Spitzer. It turns out 
that half of their sources have X-ray counterparts 
(TablelSl). Note that [BHB2007] 1 and 2 are far off- 
axis in the XMM-JVewton data and detected as a 
relatively bright unresolved X-ray source. As cor- 
roborated by near-infrared spectroscopy (Covey et 
al. 2010, subm.), sources [BHB2007] 5 and 17 
probably are background giants and remain un- 
detected in X-rays. The sources [BHB2007] 3, 4, 
8, 10, 12, 14, and 19 remain undetected in X-rays 
as well. While [BHB2007] 3, 4, 12, and 19 are lo- 
cated in chip gaps of the EPIC-PN detector, they 
are covered by the two MOS imaging detectors. 
While the sources may be hidden by higher fore- 
ground extinction, they could also be the lowest- 
mass members with the lowest X-ray luminosities. 
The fact that we detect one of the most deeply 
embedded sources, [BIIB2007] 11, indicates that 
high extinction may not be the dominant reason 
for these non-detections. 

As listed in Table |3l nine sources have suffi- 
cient X-ray counts on the EPIC-pn chip to war- 



rant a closer look at their X-ray spectra. Fits to 
these spectra were obtained using absorbed single- 
temperature APEC hot-plasma emission in Sherpa 
with a fixed setting of the metal abundances to 0.3 
times the solar value. The spectra can be reason- 
ably well fitted with such a model setup and the 
derived plasma temperatures are on the order of 
several 10^ K. The derived absorbing column den- 
sity can be empirically translated into magnitudes 
of visual extinction (see Section [OT) . The result- 
ing values indicate extinction at levels consistently 
and considerably lower than the values derived 
from the extinction map (see Table [3]). However, 
while the extinction map traces the entire extinc- 
tion along the line of sight toward the background 
sources, the X-ray extinction measures the extinc- 
tion toward the respective source. This indicates 
that the detected sources are located toward the 
front of the cloud and the undetected sources may 
indeed be hidden by higher foreground extinction. 
We have also extracted X-ray light curves for these 
sources. [BHB2007] 9, 13, and 16 - all located in 
the inner cluster - show conspicuous flares, likely 
of coronal origin, corroborating their classification 
as YSOs. 

We note that two of the brightest X-ray sources 
in the B 59 field, with more than 3000 counts each, 
are located in the off-core area. These are 2XMMi 
J171149.9-272050 and 2XMMi J171214.1-272140. 
Both have optical, near- and mid-infrared (IRAC) 
counterparts, and the latter source has a 24 /im 
counterpart. Their SEDs look like stellar pho- 
tospheres (unextincted for 2XMMi J171149.9 and 
Av 1.2 mag in the case of 2XMMi J171214.1), 
even though the 24 /xm point of XMMi J171214.1 
falls below a fit. Compared to extinction map es- 
timates Ay ~ 4.7 mag and Ay ^ 3.7 mag, the 
low extinction as derived from the SEDs indicates 
that the sources are located in front of the cloud. 
The effective temperatures derived from the SED 
fits are on the order of 3500 K, i.e., these sources 
might be active M dwarfs. 2XMMi J171149.9 
shows two flares in its X-ray lightcurve, indica- 
tive of coronal activity, while 2XMMi J171214.1 
only shows slow variation. The X-ray spectra of 
both sources indicate very low levels of extinction, 
consistent with the SED results. The X-ray spec- 
trum of 2XMMi J171214.1 can be weh fitted by 
absorbed single-temperature hot-plasma emission 
while more than one component would be needed 
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for 2XMMi J171149.9. Note that in both cases, 
also the X-ray to optical flux ratios would be con- 
sistent with an identification as M dwarfj^. A 
third bright X-ray source in the southern part of 
the B 59 field, 2XMMi J171153. 2-273607 does not 
have an optical or 2MASS counterpart and is only 
covered in two IRAC bands (1 and 3); it is also 
not covered at 24 /im (note that it seems to have 
a counterpart in the two IRAC bands). 

4. Summary 

We have used X-ray data to probe the Pipe 
Nebula complex for signs of a population of 
class III pre-main sequence stars. Such a pop- 
ulation would have been missed by our previous 
work using mid-infrared excess emission to iden- 
tify YSOs. To survey the entire complex, the only 
complete dataset is the relatively shallow ROSAT 
All-Sky Survey. Due to the nearby location of the 
Pipe Nebula, however, the absolute sensitivity of 
these data is still meaningful since we would ex- 
pect about half of a population of G and K PMS 
stars to be detected. We find a deficit rather than 
an excess of ROSAT X-ray sources in the Pipe 
Nebula relative to its immediate surroundings. 
This finding is consistent with the X-ray popula- 
tion being dominated by background sources in 
the Galactic bulge, partly extinguished by the 
foreground Pipe Nebula. Only a single ROSAT 
source has an error circle that overlaps with one 
of the previously identified extinction cores, i.e., 
at most one ROSAT source is associated to the 
on-core regions. 

In a more detailed look, we analyze three 
archival XMM-Ncwton pointings of prominent re- 
gions in the Pipe Nebula, encompassing B 59, the 
PMR, and B68. We analyze the infrared SEDs 
of all X-ray sources with near-infrared counter- 
parts {Ks < 12.4 mag) in the on-core regions 
(plus most in the off-core regions, limited by data 
availability) and find that all sources in the on- 
core regions and all but one source in the off-core 
regions can be explained by extincted stellar pho- 
tospheres. Among these three observations, the 
known cluster of YSOs in B 59 is the only re- 
gion where the density of X-ray sources in the 

^^see, e.g., lAeueros et"all ll200g) for stellar X-ray to optical 
flux ratios; with optical photom etry from the Guide Star 
Catalog 2.3 l lLasker et al.ll200Sf) 



on-core regions is indeed higher than the back- 
ground source density. Without the known YSOs, 
B 59 also shows an underdensity of sources in the 
on-core region, indicating that there is no signifi- 
cant YSO population beyond the sources that are 
already known. Roughly half of the previously 
known Spitzer-identified YSOs have X-ray coun- 
terparts with typical spectra; some show clear 
signs of variability. 

In the on-core regions of the XMM-JVewton ob- 
servations, we have found an upper limit of 6 
evolved YSOs since this is the number of X-ray 
sources with Ks-hand counterparts, not counting 
the known YSOs. However, the fact that the two 
regions with any detections in the on-core regions, 
B 59 and the PMR, both show underdensities with 
respect to the surrounding background source den- 
sities (not taking into account the known YSOs in 
B 59) suggests that a more realistic upper limit 
is closer to zero. A more stringent upper limit 
results when additionally requiring detections at 
24 /im. Only two sources in the on-core regions of 
B 59 and only four more sources in the B 59 field 
fulfill that criterion. We therefore conclude from 
the available X-ray data that there is no indica- 
tion of an extended population of class III objects 
in the Pipe Nebula, corroborating the previously 
stated lack of star formation activity. 

Facilities: XMM, Spitzer, ROSAT, 2MASS 
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Table 1 

Source selection for the three XMM-JVewtou pointings 





B59 


PMR 


B 68 


total 


Av(toty [mag] 


Av = 5.4 ± 3.5 


Av = 8.7 ± 2.6 


Av = 1.2 ± 1.2 




no. 2MASS in FOV 


10320 


19594 


13790 


43704 


no. 2MASS (Ks < 12.4 mag) 


1318 


4625 


1687 


7630 


no. 2XMMi 


101 


38 


112 


251 


no. 2XMMi w/Ks < 12.4 mag 


23 


9 


14b 


46 


no. 2XMMi w/Ks < 12.4 mag on-core 


12 


1 





13 


prev. known YSOs^ 


< 15 








< 15 


chance alignments'^ 


3.5 ± 1.9 


6.6 ± 2.6 


6.9 ± 2.6 


17.0 


chance alignments on-corc'^ 


1.7 ± 1.3 


2.6 ± 1.6 


1.0 ± 1.0 


5.3 



^Mcan Ay and standard deviation in the three fields, from extinetion map; see text. 
"^We exclude nine low-extinction (and off-core) objects of these 14 sources; sec text, 
'^as listed bv lForbrich et al.| (I2009D 

'^The chance alignments are the 2MASS sources {Kg < 12.4 mag) that would randomly fall onto the 
combined area covered by the 3cr 2XMMi search radii. 



Table 2 

2XMMl SOURCES WITH 2MASS-X5 COUNTERPARTS {K < 12.4 MAG) IN THE ON-CORE REGIONS 



No 


Field 


2XMMi 


Avltot)" 
[mag] 


Av(sub)'' 
[mag] 


srcML'' 


r(2MASS)'^ 
["] 


2MASS-Ks 
[mag] 


err 
[mag] 


ID 


1 


B59 


J171059. 9-272543 


7.3 


3.6 


9.60E+01 


2.70 


10.30 


0.02 




2 


B59 


J171103. 9-272256 


7.9 


4.2 


2.86E+03 


1.40 


7.76 


0.02 


[BHB2007] 1/2 


3 


B59 


J171113. 1-272936 


7.2 


2.6 


1.13E+01 


2.70 


9.07 


0.02 




4 


B59 


J171116. 4-272514 


15.0 


7.8 


9.77E+03 


1.41 


8.75 


0.03 


[BHB2007] 6 


5 


B59 


J171117. 5-272535 


13.0 


8.5 


2.29E+02 


1.41 


10.72 


0.02 




6 


B59 


J171121. 5-272741 


15.7 


12.4 


1.06E+03 


0.55 


8.98 


0.03 


[BHB2007] 9 


7 


B59 


J171127. 0-272348 


18.1 


15.8 


5.33E+03 


0.10 


9.08 


0.03 


[BHB2007] 13 


8 


B59 


J171129. 4-272536 


15.6 


12.4 


6.19E+02 


0.05 


10.70 


0.02 


[BHB2007] 15 


9 


B59 


J171130. 3-272629 


14.4 


10.4 


1.32E+04 


0.26 


8.89 


0.02 


[BHB2007] 16 


10 


B59 


J171141. 8-272547 


8.7 


4.5 


9.48E+03 


0.13 


8.99 


0.04 


[BHB2007] 18 


11 


B59 


J171148. 0-272633 


8.5 


4.9 


3.70E+02 


0.91 


10.87 


0.03 




12 


B59 


J171231. 5-271847 


5.3 


2.5 


1.76E+01 


2.06 


11.64 


0.03 




13 


Ring 


J173402. 6-255219 


10.2 


3.0 


8.67E+00 


1.88 


12.21 


0.03 




^visual extinction from extinction maps with the total extinction, Ay 
lAlves et al.ll2007D 

''2XMMi detection likehhood. seelWatson et al.l 1200911 


(tot), and with large-scale structure removed, Av(sub) 



"^distance to nearest 2MASS source 
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Table 3 

x-ray properties of ysos in b 59 



[BHB2007] total counts kT Nh red. Av(X)" Av^(lot) 

[keV] 10^^ /cm^ [mag] 



1&2 


662'' 




3 


nd 




4 


nd 




5 


nd 




6 


2032 


1. 


7 

8 


343 
nd" 


8. 


9 


391 


3. 


10 


nd 




11 






12 


nd 




13 


1289 


1. 


14 


nd 




15 


285 


0. 


16 


2700 


1. 


17 


nd 




18 


2211 


1. 


19 


nd 




20 


114'' 





mae 



0-63l|l;o3 0-^2 3.2±0.2 15.0 
■ 97^2^1 1.18+° 0.56 6.0±0.7 15.0 

^to'.ll l-93j:o:i6 0-45 9.7±0.9 15.7 



4+o;°7 0.67j:°;°3 0.88 3.4±0.2 18.1 



0.86Inn„ 0.40 4.3±0.3 15.6 



O + 0.09 „ o(. + 0.07 

^ — 0.09 — 0.06 

Slo.os 0.53j:°;°^ 0.62 2.7±0.1 14.4 

2lo.o2 0.04j:°°J 1.09 1.4±0.4 8.7 



^optical extinction from the fit to the X-ray spectrum, Av(X); and from the total- 
extinction map, Av(tot). 

''blend, detected in EPIC-MOS only due to large off-axis angle; best fitted with two 

separate, extinctcd plasma components 

'^possibly a weak X-ray source, but not picked up in source detection 
''poor fit; few data points 
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3:00:00.0 2:00:00.0 1:00:00.0 0:00*0.0 359:00:00.0 358:00:00.0 357:00:00.0 



Galactic longitude 



Fig. 1. — A background-subtracted extinction map (jAlves et al.l 120071 ) of the Pipe Nebula with ROSAT 
Ah-Sky Survey sources shown as small circles. Three large circles mark the areas studied in more detail 
using XMM-ATewto n indicating the usab le fields of view. The continuous line delineates the area surveyed 
with Spitzer-MlPS (jPorbrich et al.ll2009[ ). 
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Fig. 2. — ROSAT hardness diagram (HR l and HR2) for the sources in the Pipe Nebula region, 'x' with 
sohd hnes, and Taurus PMS sources from iKenvon et aL ( 2008h . dashed hnes. Data only shown for sources 
that have both hardness ratios, HRl and HR2, with errors < 0.3. 
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Fig. 3. — XMM-Newton view of B59, using merged data of the PN and MOSl/2 detectors. Contour lines 

show the Av = 1.2 mag levels from the background-subtracted extinction map, delineating the "on-core" 
regions. Circles show 2XMMi sources with 2MASS counterparts having iCs-band detections (< 12.4 mag). 
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30.0 17:35:00.0 30.0 34:00.0 33:30.0 

Right ascension 

Fig. 4. — ^MM-Newton views of the Pipe molecular ring, using merged data of the PN and MOSl/2 
detectors. Contour hues and symbols as in Fig. [3l 
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Fig. 5. — XMM-Newton views of B68, using merged data of the PN and MOSl/2 detectors. Contour 
lines and symbols as in Fig. |3l Note that only the final selection of sources is shown. Among the three 
observations, this field has by far the lowest extinction in the off-core regions (see Table [1]). 
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